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The work of the Lubhn group on the non-equilhbrium transport through the quantum dot coupled 
to external leads (normal or superconducting) and subject to external time dependent fields has been 
reviewed. 



I. INTRODUCTION 

This paper is an extended version of the talk given at the LFPPI-network conference in Lodz 
(Poland). It makes extensive use of our published work In particular all figures which appear 

in the present review have already been presented there. 

Due to recent developments in fabrication of small electronic devices equilibrium and non- 
equilibrium properties of mesoscopic systems have been of great interest during last decade. In 
particular a lot of work has been done on electron transport through quantum dots. 

The quantum dots are small pieces of matter containing few to few tens of electrons. Very often 
they are electrically defined in a two dimensional electron gas placed at the interface in GaAs-GaAlAs 
heterostructure. The important point about quantum dots is that they can be viewed like impurities 
in solids but at the same time their coupling to external environment can be controlled and the 
noneguilibrium transport can be studied. 

The discovery of the Kondo effect in quantum dots [l^, ll^ connected to external leads by tun- 
nel junctions has resulted in increased experimental and theoretical interest. The Kondo effect in 
quantum dots manifests itself as an increased conductance G of the system at temperatures lower 
than the Kondo temperature Tk- It is due to formation of the so-called Abrikosov-Suhl or Kondo 
resonance at the Fermi energy. This is a many-body singlet state involving spin on the quantum dot 
and electrons in the leads. 

Here we review our studies of stationary nonequilibrium transport with help of Keldysh Green's 
function formalism. To obtain the information about the time dependent phenomena in a quantum 
dots which parameters do depend on time we calculate the time evolution operator and then time 
dependent charge and current. 

II. TIME-INDEPENDENT TRANSPORT 

In this section we shall present the results for time independent transport through quantum dot 
with a single energy level. We start with asymmetrically coupled quantum dot (section II. A), discuss 
the role of Van Hove singularity in the density of states of the electrodes (section II. B) and consider 
the dot coupled to one normal and other superconducting lead (section II. C). 



2 




0.00 I ' ' ' ' ' ' I 

-4-3-2-10123 



(0 



FIG. 1: The equilibrium density of states on the quantum dot obtained within EOM (solid line) and NCA (dashed line). 
Note the relative shift of the spectral weight with respect to the chemical potential ji = Q which results in different occupations 
shown in the inset. The parameters are: Ed ~ —2, Vi^ — F^j = 1 and T = 10~^. 



A. Non-equilibrium Kondo effect in asymmetrically coupled quantum dots 

It is well known that in systems containing quantum dot the Kondo effect manifests itself (at low 
temperatures) as an enhancement of the conductance at zero source-drain voltage[10], Vsd = 0. 
Occasionally the Kondo peak in conductance appearing at non-zero voltages Vsd 7^ ^ has been 
observed 0, 0|, the so called anomalous Kondo effect. This phenomenon has been investigated 
experimentally for the dot coupled weakly to one and strongly to the other lead. It was observed 
that the source-drain voltage, at which the peak appeared, scales roughly linearly with a gate voltage 
Vg. Performing model calculations based on non-equilibrium transport theory we have found that 
the emergence of the Kondo peak at non-zero voltages is caused by asymmetric coupling to external 
leads. 
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FIG. 2: The non-equilibrium density of states obtained within a) - EOM and b) - NCA for the symmetric Tl ~ Fjj (solid 
lines) and asymmetrically coupled quantum dot with F_l = 2Fj{ (dashed) and Fl = ^F_r (dotted lines), — —/il ~ 0.2 and 
the other parameters are the same as in Fig. 



The system under consideration is described by single impurity Anderson model and the Landauer- 
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type formula for the current flowing in the system, left lead - quantum dot - right lead is 11 
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where fL/ni^) denotes the Fermi distribution function, G^^{e) - retarded on-dot Green's function and 
e) - effective coupling of the localized electrons to the conduction band. 
The retarded Green's function required for the calculations of the tunneling current and the differ- 
ential conductance of the system is obtained within the equation of motion [EOM) technique with 
slave boson representation of the electron operators and non-crossing approximation (NCA). In 
both cases we assume U —>■ oo. As we have checked, the EOM method gives correct position of the 
Kondo peak, however it leads to incorrect width of the peak and occupation on the dot. Therefore 
we have used NCA technique which is generally accepted for description of the system in the Kondo 
regime. 

In the Fig^we present the equilibrium density of states {DOS) of the quantum dot coupled to 
two leads obtained by means of the NCA and EOM approaches. The main features of the DOS 
remain the same in both cases. However the height and width of the Kondo peak are much larger 
in NCA. 

FigI21gives the non-equilibrium DOS obtained for symmetric {Tl = Tr) and asymmetric couplings 
and Tl = ^Tj 
those of the Fermi levels of the leads 



Tl = 2Tji and Tl = ^Tj?- We see that the Kondo peak is always located at energies coinciding with 
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Thus in non-equilibrium we get two Kondo resonances in 
the DOS pinned to the Fermi energies of corresponding leads. 
The differential conductance spectrum of the same system is shown in the FigOl For comparison 
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FIG. 3: The difFerential conductance {G(eVsD) ~ dJ/d{eVsD)) obtained within a) - EOM and b) - NCA for the symmetric 
Tl ~ Tj; (sohd hnes) and asymmetrically coupled quantum dot with Fl = 2r_R (dashed) and Fl ~ ^Fr (dotted lines). 0] 



we have also plotted the conductance of the symmetrically coupled quantum dot. In the symmetric 
case (Tl = Tfi) the Kondo peak is located exactly at zero bias (Vsd = 0) but for > F^j (F^ < Tr) 
it is shifted to the negative (positive) voltages. This finding is in nice qualitative agreement with 
experimental situation in transport through the quantum dot in the presence of asymmetric barriers. 
Further details regarding this anomalous Kondo effect can be found in |2i. 



B. Van Hove singularities and the Kondo effect 

Here we present some results regarding influence of the band structure of the leads, in particular 
Van Hove singularities in the DOS on the tunneling current through the quantum dot. We have 
assumed the conduction electron energies of the leads in the form characteristic for two-dimensional 
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tight-binding spectrum which is known to possess the logarithmic singularity in the middle of the 
band. To calculate the retarded Green's function we have used Hubbard I type decoupling and 
performed calculations of the density of states (DOS) on the dot and the differential conductunce 
trough it. To see the influence of spectrum of the electrode on transport we repeated the calculations 
for constant, Lorentzian and Van Hove singular DOS. 
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FIG. 4: The density of states on the dot calculated for V — O.IW, T — 10~'*Vy (W is the bandwidth) and for three different 
spectra of electrons in the leads. Note the structure due to Van Elove singularities and robustness of the coherent (Kondo) 
feature at the Fermi level. 

Fig. m shows the effect of the Van Hove singularity in the leads on the density of states of the 
electrons on the dot and figure shows the corresponding differential conductance spectra. One 
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FIG. 5: The differential conductance G (in units of e'^/h) for three spectra of electrons in the leads. The Van Hove singularity 
affects the conductance and produces additional structure around the Fermi level. 



observes marked changes of the tunneling current and the differential conductance due to the Van 
Hove singularities (now shifted by ±0.0514^). The dip seen below Fermi level {fiL — /Wi? < 0) is due 
to the real part of the self-energy, while the Kondo like structure above the Fermi level (solid curve) 
is connected with the imaginary part, which is proportional to the lead DOS. The structure thus 
reflects the presence of Van Hove singularity in the spectrum of carriers in the leads. 
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C. Electron transport through a quantum dot coupled to normal metal and superconductor 

Here we present some of the results of our studies concerning transport properties of the quantum 
dot coupled to the normal metal and i^CS"- like superconductor (A^ — QD — S) in the presence of the 
Kondo effect and Andreev scattering. The system is described by single impurity Anderson model 
in the limit of strong on-dot interaction. We have used EOM method to calculate non-equilibrium 
Green's functions together with the modified slave boson technique. The details of calculations can 
be found in jl] and here we present the results for the quantum dot DOS only, as it gives valuable 
information about the system. 

In Fig. iniwe show the DOS of the quantum dot for various positions of the dot energy level Ed- 
Clearly, the Kondo effect manifests itself in the resonance at the Fermi energy which survives the 
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FIG. 6: The density of states of the quantum dot in N ~ QD — 5* system for various values of the dot energy level Ed- Other 
parameters are following: Tjv = Tg = 0.02, A = 0.1, eV = 0, T = 10~^ in the units of the bandwidth W. ||| 

presence of superconductivity in one electrode. The additional structure at a; = ±A = ±0.1 reflects 
the superconducting gap A. 

In the next step we consider DOS in the non-equilibrium situation {eV = fij\f — Recall 
that in the A^ — QD — N system for eV ^ there emerge two resonances pinned to Fermi levels 
of the left and right leads, respectively. In the present case there is a gap in the spectrum of the 
superconducting lead and we observe only one resonance pinned to the normal metal Fermi level 
(/iAr) (see Fig. [7j). The presence of one superconducting electrode opens a possibility for the Andreev 
reflections. This is process in which an electron coming to the N-S interface from the normal metal 
side is reflected back into metal as hole, while the Cooper pair enters the superconductor. From the 
general expression for the nonequilibrium current (P) in the A^ — QD — S system one expects the 
Kondo peak in the tunneling mediated by the Andreev reflections. The corresponding part of the 
total current is written in the form Q 

Ja = r ^T^si^)[f{u^ - eV) - fiu; + eV)], (2) 

where we have introduced "transmittance" T^g{uj) associated with the Andreev tunneling. At 
zero temperature and at energies less than the superconducting gap, T^g{uj) can be regarded as a 
total transmittance because Andreev tunneling is the only process allowed in these circumstances. 
T^g{u) for different values of the eV is plotted in the Fig|Hl The broad resonances at a; ~ ±0.06 are 
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FIG. 7: Equilibrium {eV = 0) and non-equilibrium (eV = ±0.01) density of states of the quantum dot in the A'^ — QD — S 
system. Other parameters have following values: Tn = Fs = 0.02, A = 0.1, Ed = —0.08, T = in units of the bandwidth 
W. i 
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FIG. 8: TjYg(a') for different values of the bias voltage eV — (solid line), —0.03 (dashed) and 0.03 (dotted line). Ed 
Fiv = Fs = 0.01 and A = 0.1. Inset; large scale view of the equilibrium T^g^ui) Q]. 
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reflections of the dot energy level = 0.08 for electrons and holes, shifted from its original position 
due to renormalization caused by the strong Coulomb interaction. The important point is that there 
is no Kondo peak in equilibrium {eV = 0) transmittance. However as soon as we go away from the 
eV = 0, we can observe the Kondo peaks at energies u = ±eV and the strong asymmetry between 
negative (dashed line) and positive (dotted line) voltages. While in the former case we have very 
well resolved resonances, in the later these resonances are strongly suppressed. This asymmetry is 
strictly related to the density of states (see Fig. ©), where we also observe such asymmetry. 

Since equilibrium transmittance T^g^u) does not show the Kondo peak, we cannot expect it in 
the differential conductance GAieVso) = dJAl d^eVsn) with J a deflned by It turns out that 
GA{,eVsD) is very sensitive to the position of the dot energy level (see Fig. 7 in ref. @|). More detailed 
discussion of these results can be found in Ref. 
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III. TIME-DEPENDENT TRANSPORT 



In this section we consider the transport properties of a quantum dot (QD) under the influence of 
external time-dependent fields. The high-frequency signals may be applied to a QD and the time- 
dependent fields will modify the tunneling current. New effects have been observed and theoretically 
described, e.g. photon-assisted tunneling through small quan tum dots with well-resolved discrete 
energy states 



LOton-assistea tunnenng tnrougn small quan tum aots witn weii-resoivea aiscrete 
17, 3], photon-electron pumps 2l| and others. One can investigate the 



current flowing through a QD under periodic modulation of the QD electronic structure j22] or 
periodic (non-periodic) modulation of the tunneling barriers ^3] and electron energy levels in both 
(left and right) electron reservoirs [23^ 

We consider the simplest case of the quantum dot with multiple energy levels without the electron- 
electron Coulomb interaction coupled to two leads (right and left) in the presence of external mi- 
crowave (MW) fields which act on all parts of the system. Under the adiabatic approximation 
the single-electron energies of the time dependent driven system can be represented in the form 
£j:{t) = + A cos cut. We describe the dynamical evolution of the charge localized on the QD and 
the current flowing through the system in terms of the time evolution operator (for details see Refs 
[1,3,9,24]. 

In the first subsection we consider the influence of the band structure of the leads on electron 
transport through the QD and in the next subsection we consider the time-dependent transport 
through the QD with additional tunneling channel. 



A. Band structure effects in time-dependent electron transport through the QD 

Here we are going to study the influence of the singularities of the leads band structure and 
therefore we are forced to go beyond WBL. Within the formalism of the non-equilibrium Green's 
functions it is extremely difficult to calculate, e.g. the time dependent current jL{t) tunneling from 
the left lead to the QD. We recall that the general formula for the current flowing e.g. from the 
left lead into the QD is usually expressed in terms of the QD retarded Green function G^(t,ti). 
However, G^(t, ti) satisfies the Dyson equation with subsequent double time integrations and for the 
time-dependent Hamiltonian it is a very difficult task to calculate it. Therefore, in the following, 
we calculate the QD charge and current using the evolution operator method. The potential drop 
between the left and right leads is given by /il — /i/? = eVs-d and Vs-d is the measured voltage 
between source and drain. In experiments the gate voltage controls the position of QD's energy level 
ed and to mimic measurements of the QD charge or current vs gate voltage we have calculated them 
against e^. To study the effect of the leads band structure on the electron transport through the 
QD, we have assumed the two-dimensional tight binding simple cubic crystal spectrum (2D-TB) for 
leads conduction electrons, which is known to possess the logarithmic singularity in the middle of 
the band. 

In Fig. ini in the left panels we show the time-averaged QD charge, tunneling current and the 
derivative of the tunneling current with respect to the QD energy level e^. In the right panels the 
corresponding differences {A{n{t)), A{jL(t)) and A.{dji) / ded) between the results obtained for 2D- 
TB leads DOS and obtained within wide-band limit (WBL) usually used in literature are given. 
The dashed lines correspond to the static case [u = 0) and the solid lines correspond to the time- 
dependent transport. 

One can see, that, e.g. for = 1, the tunneling current calculated for 2D-TB leads DOS is greater 
than that obtained within WBL up to a factor of ~ 0.30 whilst for = 1 up to ~ 1.50 for the 
static case and this difference decreases with increasing amplitude (A^ = Al/2, Ar = 0). It is 
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FIG. 9: The time averaged QD charge (n(i)), tunnehng current {jLit)) and derivative d{jL{t)} /ded (in arb. units) vs ed obtained 
for 2D-TB leads DOS - the left panels. In the right panels the corresponding differences between the results showed in left 
panels and those obtained within WBL {A{n{t)), A{jL{t)) and A{djL) /ded) are given. The thick (thin) curves correspond to 
Al = 8, Ad — 4, Ah = (Al ~ 4 ,A£j — 2, A^ — 0) and uj = 2, hr — 0, fiL = 4. The dashed curves correspond to the 
time-independent case. |^ 

interesting that the A{jL{t)) curves for different exhibit local minima or maxima at the same 
ed- For example, A{jL(t)) for the static case possesses a local minimum for = 2 (the middle 
point between the chemical potentials fii and fin) but for the time- dependent case we observe local 
maxima, greater for Al = 4: and lesser for greater amplitude = 8. The bottom panels show the 
derivative d{jL{t)) /ded vs ed and the corresponding differences obtained for two considered models 
of the leads DOS. As before, the greatest differences are observed for the static case (dashed lines) 
for QD energy levels equal to the chemical potentials fi^ or fin. For the time-dependent transport 
these differences are much smaller and hardly depend on amplitudes A^ and A^, contrary to A(ji(t)) 
which are greater for smaller values of Ai, and A^. 

In the next step we consider the influence of the lead DOS structure on the photon-assisted 
tunneling through a quantum dot. In experiment, in average current vs gate voltage curve, a 
"shoulder" is observed on the left side of the main resonant peak for the case of a MW field applied 
on one lead only. In Fig^J we compare the structure of the main resonant peak in the average 
current {iiif)) and the derivative d{ji,)/ded plotted vs obtained within WBL (right panels) and 
calculated for 2D-TB leads DOS (left panels) at half-filled bands. The remarkable differences are 
visible for negative values of e^. The main resonant peaks are broader and lower and their height 
is almost independent on the amplitude A^. The photon-assisted tunneling for = ±1 is clearly 
visible (lower panels of FigHUJ. 
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FIG. 10: Left panels - the average current {jhit)} and derivative d{jL{t)} /dtd vs ed for external field applied only to the right 
lead in the case of the 2D-TB leads DOS. The solid curves correspond to Ar/uj — 0.3 (thin lines), 0.5 and 0.7 (thick lines). 
The dotted curve is the resonant peak without external field. The parameters are Fl ~ Tr = O.lcu and hl ~ —ftR = 0.05aj. 
For comparison, on the right panels the same as on the left but obtained within WBL. |^ 



B. Electron transport through a QD in the presence of a direct tunneling between leads 

The progress of nanomaterials science has enabled the experimental study of the phase coherence of 
the charge carriers in many mesoscopic systems. Especially, the Fano resonances in the conductance 
were observed [1^ , which imply that there are two paths for transfer of electrons between a source and 
a drain. The recent experimental and theoretical study with a low-temperature scanning tunneling 
microscope (STM) of the single magnetic atom deposited on a metallic surface showed also the 
asymmetric Fano resonances in the tunneling spectra jl^ 22, H^. The STM measurements 
indicate that in tunneling of electrons between STM tip and a surface with a single impurity atom 
two different paths are present. The electrons can tunnel between the tip and the adsorbate state 
and directly between the tip and the metal surface. 

In this subsection we address the issue of a QD with a bridge channel between a source and a 
drain driven out of equilibrium by means of a dc voltage bias and additional time-dependent external 
fields. 

We model the QD coupled to the left and right electron reservoirs with the additional bridge 
tunneling channel between them (represented by the coupling Vl/j) by the usually used Anderson- 
type Hamiltonian where U = 0. 

Here we show the numerical results of the time-averaged current {jL{t)) and its derivatives with 
respect to the QD energy level position and the chemical potential fiL (or equivalently, with respect 
to the gate and source-drain voltages) for different sets of parameters which characterize our system. 

At first, we consider the dependence of {jhif)) vs. Sd for given values of the left lead chemical 
potential /i^. In Fig. [TTlwe present such curves for different values of hl obtained for hl = 8. The 
left (right) panel corresponds to Vrl = (Vrx = 10). In the case of vanishing over-dot tunneling 
channel (the left panel) the current has a simple structure - a single peak localized in the middle 
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FIG. 11: The average current {jL{t)) against Ed for given values of hl (beginning from jiL = —4 up to hl = 8). The left and 
right panels correspond to Vrl = and Vrl = 10, respectively, fj,R = 0, V — 4, Al ~ 2, A^j — 1, Ar = 0, uj = 2. 
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FIG. 12: The derivatives of the average current against Sd with respect to the QD energy level Ed, d{jL{t)) /dsd, for given values 
of fiL (beginning from fiL = —4 up to fiL = 8). The broken curves correspond to hl ~ —4. The left (right) panel corresponds 
to Vrl = {Vrl = 10) and the other parameters as in Fig. 0. 

between fi^ and fiR for smaller values of /i^. The width of this peak increases with increasing and 
for greater values of the current is almost independent of Ed localized inside the energy region 
between fi^ and fii- For non- vanishing over-dot tunneling (Fig. ^2 the right panel) the curves 
become asymmetric. With increasing source-drain bias, the current possesses greater values 
in comparison with the Vrl = case due to the direct tunneling between both leads. Note, however, 
that due to the interference effects the resulting (jL(t)) curves are asymmetric. The interference 
effects are most visible for Ed lying approximately in the region {fiR,fiL)- 

In Fig. ^Jwe show the derivatives of the average current vs. the QD energy level Ed obtained for 
some values of fii. There are the results of the differentiation of curves shown in Fig. ^2 Again, 
the most visible differences between the results obtained for Vrl = and Vrl 7^ are present for 
the QD energy level Ed localized approximately between chemical potentials of both leads (compare, 
for example, the curves calculated for /x^ = 8). 

Fig. Unipresents the comparison of the d{jL(t)) /d^L vs. fiL curves calculated for vanishing Vrl (left 
panels) and for Vrl = 10 (right panels) for two different values of the amplitudes Al (A^ = Al/2, 
= 0). At the vanishing value of Vrl, the shape of the curves is symmetrical in relation to the 
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FIG. 13: The derivatives of the average current against fiL with respect to fiL, d{jL{t)) /dfiL, for given values of Sd (beginning 
from Ed ~ —4 up to Ed = 8). The left (right) panels correspond to Vrl = {Vrl ~ 10) and upper (lower) panels correspond 
to Al = 2 (Al = 4). The other parameters as in Fig. ^ 

values fiL = £d although for greater some shoulders appear on both sides of the corresponding 
peaks in the distance ~ from the curve centers. For non- vanishing Vrl, the corresponding 

curves are approximately asymmetric and for large values of fii they tend to constant, non-zero 
values corresponding to linear increasing of the current at large /x/,. It is interesting that with the 
increasing amplitudes Al and very clear structures appear on both sides of the corresponding 
curves. Note, that all these curves can be obtained, for example, from the one calculated for £d = 
and moved along the /iL-axis by the corresponding value (equal to Ed)- 

Fig. El presents the average current {jL{t)) obtained for the case in which only the QD energy 
level Ed oscillates with some frequency {u > T ) and for the small source-drain voltage Hl — I^r = 0.2. 
For vanishing V/ji (broken curves) we observe for small amplitude A^ only a central resonant peak 
(Fig. m^). With increasing A^, the subsequent peaks appear and the distance between them and 
the central peak is an integer multiple of the frequency uj (side bands). The location of peaks is 
independent of the amplitude A^ but their relative intensity values change and with increasing A^ 
the height of the central peak is reduced. If we take the additional over-dot tunneling channel into 
consideration, especially for small A^, the asymmetric shape of the current curve is observed and 
this asymmetry increases with the increasing strength of the over-dot coupling between both leads. 
With the increasing amplitude A^; this asymmetry is reduced largely due to the extra, photon-assisted 
tunneling peaks whose strength increases with the increasing A^. 

As a last problem we consider the QD with the periodic rectangular-pulse external field applied to 
each QD-lead barrier. We assume that the influence of these fields on the system under consideration 
is equivalent to the following time-dependence of the matrix elements Vl and Vr corresponding to 
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FIG. 14: The average current {jhit)) against Sd for the oscillating QD energy level at Vrl = 0, 4, 7, 10 - broken, thin, thick 
and very thick curves, respectively. The panels a, b, c and d correspond to Ad = 1,3,5 and 7, respectively, a; = 5, F = 1, 
V = 4, /ii = 0.2, Al = Afl = 0. II 

the coupling of the QD with both leads: 

VL{t) = ^,VR = Qfor{]<t<T/2 
VL{t)=Q,VR = A for T/2<t<T. 

The bridge tunneling channel is present all the time and the calculations were performed for 
Vrl = 4 and Vrl = 10. In the upper part of Fig. UHlwe show the time- dependence of VLit). During 
the periodic coupling of the QD to the left or right lead the QD charge oscillates and the amplitude 
of these oscillations decreases with the increasing coupling Vrl between both leads. Note, that the 
greater Vrl enhances the equilibrium value of the QD charge at Vl = and 7^ (in the first 
half-period), but reduces it for Vl 7^ and Vr = (i.e. in the second half-period). In two lower parts 
of Fig. El we show the current leaving the right and left leads for two values of the coupling Vrl- 
The effect of switching-on and -out of the QD-leads coupling on the currents is markedly visible. 
The time after which the currents achieve their steady values hardly depends on the coupling Vrl, 
although the width of the corresponding current peaks is smaller for lower values of Vrl- 
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FIG. 15: The time-dependence of n{t), I hit) and lR{t) for the case when the QD-lead barriers are changing in time as shown 
for Vhit) in the upper panel (l^fl(t) is out of phase with a phase difference of tt) for two values of the over-dot bridge tunneling 
channel strength corresponding to Vrl ~ 4 (thin lines) and Vrl ~ 10 (thick lines), sa ~ 1, V = A, hl = i, Hr = 0. Q 

IV. CONCLUDING REMARKS: IS ALL THIS RELATED TO QUANTUM INFORMATION 

PROCESSING? 



The problems discussed here are of importance for small structures working at low temperatures. 
In such cases the presence of even single electron on the structure prevents other electrons to tunnel 
onto it. This is because small capacitance of these systems results in large charging energies which 
have been parametrised here by the Coulomb on dot repulsion. 

As the ultimate goal of the miniaturisation is to build a quantum computer it maybe of interest 
to ask the question whether the systems considered will ever be used for this purpose. The answer 
is affirmative. It has been proposed to use charge [2^ or spin js^] on a quantum dot to build a 
qubit - an elementary quantum unit of information. The subsequent work has shown that the 
decoherence of a system consisting of two lateral quantum dots due to cotunneling does depend on 
the difference between chemical potentials of the two leads i.e. on the voltage. The future work will 
show whether the control of decoherence by the external bias will also be effective in more complex 
geometries. 



14 



Acknowledgements: This work has been partially supported by the grant no. PBZ-MIN- 
008/P03/2003 . 



[1] Kwapinski T, Taranko R., Time- dependent transport through a quantum dot with the over- dot(bridge) additional tunneling 

channel Physica E 18 (2003) 402-411 
[2] Krawiec M, Wysokiriski KL, Nonequilibnum Kondo effect in asymmetrically coupled quantum dots Phys. Rev. B 66 165408 

(2002). 

[3] Kwapinski T, Taranko R, Taranko E., Band structure effects in time- dependent electron transport through the quantum dot 

Phys. Rev. B 66 035315 (2002). 
[4] Kwapinski T, Taranko R, Taranko E., Time- dependent transport through a quantum dot: Mean-field, treatment of the 

Coulomb interaction Acta Phys. Pol. A 99 (2001) 293-305. 
[5] Krawiec M, Wysokinski KL, Charge on the quantum dot in the presence of tunneling current Solid State Commun 115 

(2000) 141-144. 

[6] Krawiec M, Wysokinski KL, Spectral functions of the quantum dot coupled to normal and/or superconducting leads Acta 

Phys. Pol. A 97 (2000) 197-200. 
[7] Krawiec M, Domariski T, Wysokinski KL, Do Van Hove singularities in leads influence tunneling current through quantum 

dot? Acta Phys. Pol. A 94 (1998) 411-414. 
[8] M. Krawiec, K. I. Wysokinski, Electron transport through strongly interacting quantum dot coupled to normal metal and 

superconductor Supercond. Sci. Technol. 17, 103 (2004). 
[9] Taranko R., Kwapinski T., Taranko E., Influence of microwave fields on the electron transport through a quantum dot in 

the presence of a direct tunneling between leads arXiv:cond-mat/0304121 (2003) unpublished. 
[10] T. K. Ng and P. A. Lee, Phys. Rev. 61, 1768 (1989); L. L Glazman, and M. E. Raikh, Pis'ma Zh. Eksp. Teor. Fiz. 48, 

378 (1988) [Engl, transl. JETP] Lett. 47, 452 (1988)]; S. Herschfield, J. H. Davies, and J. W. Wilkins, Phys. Rev. Lett. 

67, 3720 (1991). 

[11] Y. Meir N. S. Wingreen, and P. A. Lee, Phys. Rev. Lett. 70, 2601 (1993); N. S. Wingreen, and Y. Meir, Phys. Rev. B 49, 

11040 (1994). Y. Meir N. S. Wingreen, P. A. Lee, Phys. Rev. Lett. 66, 3048 (1991). 
[12] D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. Abush-Maggder, U. Meirav, M. A. Kastner, Nature 391, 156 (1998). 
[13] S. M. Cronenwett, T. H. Oosterkamp, L. P. Kouwenhoven, Science 281, 540 (1998). 
[14] J. Schmid, J. Weis, K. Eberl, K. von Klitzing, Physica B256-258, 182 (1998). 

[15] F. Simmel, R. H. Blick, J. P. Kotthaus, W. Wegsheider, M. Bichler, Phys. Rev. Lett. 83, 804, (1999). 
[16] T. H. Oosterkamp, L.P. Kouwenhoven, A.E.A. Koden, N.C. van der Vaart, C.J.P.M. Harmans, Phys. Rev. Lett. 78, 1536 
(1997). 

[17] Qing-feng Sun, Tsung-han Lin, Phys. Rev. B56 3591 (1997). 

[18] Qing-feng Sun, Jian Wang, Tsung-han Lin, Phys. Rev. B58, 13007 (1998). 

[19] C.A. Stafford, N.S. Wingreen, Phys. Rev. Lett. 76, 1916 (1996). 

[20] Qing-feng Sun, Tsung-han Lin, J. Phys.: Condens. Matter 9, 3043 (1997). 

[21] L.P. Kouwenhoven, A.T. Johnson, N.C. van der Vaart, A. van der Enden, C.J.P.M. Harmans, C.T. Foxon, Z. Phys. B, 

Condens. Matter 85, 381 (1991). 
[22] H.-K. Zhao, J. Wang, Eur. Phys. J. B9, 513 (1999). 
[23] Y. Goldin, Y. Avishai, Phys. Rev. B61, 16750 (2000). 

[24] M. Tsukada, N. Shima, in: "Dynamical Processes and Ordering on Solid Surfaces", Eds. A. Yoshimori and M. Tsukada 
(Springer, Berlin 1995), p. 34; T.B. Grimley, V.C.J. Bhasu, K.L. Sebastian, Surf. Sci. 121, 305 (1983). 

[25] J. Gores, D. Goldhaber-Gordon, S. Heemeyer, M.A. Kastner, H. Shtrikman, D. Mahalu, U. Meirav, Phys. Rev. B62, 2188 
(2000). 

[26] M. Plihal, J.W. Gadzuk, Phys. Rev. B63, 085404 (2001). 

[27] V. Madhavan, W. Chen, T. Jamneala, M.F. Crommie, N.S. Wingreen, Science 280,567 (1998). 
[28] O. Ujsaghy, J. Kroha, L. Szunyogh, A. Zawadowski, Phys. Rev. Lett. 85, 2557 (2000). 
[29] J.A. Brum, P. Hawrylak, Superlattices and Microstructures, 22 431 (1997). 
[30] D. Loss, D.P. Di Vincezo, Phys. Rev. A57 120 (1998). 

[31] U. Hartmann, F.K. W ilhelm Decoherence of charge states in double quantum dots due to cotunnelmg 
|arXiv:cond-mat/0208220| 



